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Abstract 

A better process understanding of how water erosion influences the redistribution of soil organic carbon (SOC) is sorely 
needed to unravel the role of soil erosion for the carbon (C) budget from local to global scales. The main objective of this 
study was to determine SOC redistribution and the complete C budget of a loess soil affected by water erosion. We 
measured fluxes of SOC, dissolved organic C (DOC) and CO2 in a pseudo-replicated rainfall-simulation experiment. We 
characterized different C fractions in soils and redistributed sediments using density fractionation and determined C 
enrichment ratios (CER) in the transported sediments. Erosion, transport and subsequent deposition resulted in significantly 
higher CER of the sediments exported ranging between 1.3 and 4.0. In the exported sediments, C contents (mg per g soil) of 
particulate organic C (POC, C not bound to soil minerals) and mineral-associated organic C (IVIOC) were both significantly 
higher than those of non-eroded soils indicating that water erosion resulted in losses of C-enriched material both in forms 
of POC and IVIOC. The averaged SOC fluxes as particles (4.7 g C m~^ yr~^) were 18 times larger than DOC fluxes. Cumulative 
emission of soil CO2 slightly decreased at the erosion zone while increased by 56% and 27% at the transport and 
depositional zone, respectively, in comparison to non-eroded soil. Overall, CO2 emission is the predominant form of C loss 
contributing to about 90.5% of total erosion-induced C losses in our 4-month experiment, which were equal to 18 g C m~^. 
Nevertheless, only 1.5% of the total redistributed C was mineralized to CO2 indicating a large stabilization after deposition. 
Our study also underlines the importance of C losses by particles and as DOC for understanding the effects of water erosion 
on the C balance at the interface of terrestrial and aquatic ecosystems. 



Citation: Wang X, Cammeraat ELH, Romeijn P, Kalbitz K (2014) Soil Organic Carbon Redistribution by Water Erosion - The Role of CO2 Emissions for the Carbon 
Budget. PLoS ONE 9(5): e96299. doi:10.1 371/journal.pone.0096299 

Editor: Ben Bond-Lamberty, DOE Pacific Northwest National Laboratory, United States of America 
Received October 3, 2013; Accepted April 7, 2014; Published May 6, 2014 

Copyright: © 2014 Wang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by China Scholarship Council (CSC) and University of Amsterdam (UvA). The funders had no role in study design, data 
collection and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 

* E-mail: X.Wang@uva.nl 



Introduction 

Climate change will likely modify current precipitation regimes 
influencing the global carbon (C) cycle in relation to erosion 
processes [1,2]. The length and intensity of droughts and the 
intensity of more sporadic rainfall events are predicted to increase 
for Western Europe [3], which will accelerate soil erosion. Soil 
erosion has significant impacts on the redistribution and transfor- 
mation of soil organic carbon (SOC) within a landscape [4,5]. 
Even now, there is no consensus whether soil erosion is acting as a 
net C sink [5,6] or source [7] of atmospheric COj. Therefore, 
quantitative assessments of soil organic C redistribution along 
geomorphic gradients and the processes involved become increas- 
ingly important in a changing climate to resolve this controversy 
[8] . It is crucial that such studies comprise the different processes 
associated with the redistribution of C along the slope including 
CO2 emissions as a result of changes in C mineralization upon 
erosion, transport and subsequent deposition. Based on such 
studies, complete C budgets of soils affected by erosion processes 
can be determined. 

Soil erosion seems to preferentially remove fresh and more 
labile materials from C rich topsoUs in upslope eroding positions, 
i.e. SOC with low density (e.g. free light fraction) and dissolved 
organic C (DOC) [7-10]. However, the fate of this organic C has 
rarely been studied. It is well known that most of the eroded 



sediments are re-deposited close to the source areas and in the 
catchment (e.g. [4,1 1]). Deposition of C enriched sediments lead to 
accumulation of SOC in the downslope positions. The eroded and 
deposited C can be stabilized by interaction with minerals thereby 
decreasing mineralization of deposited C in soil profdes [12]. In 
addition, soil erosion could affect dissolved organic carbon (DOC) 
dynamics in soils. Wang et al. [12] found higher DOC 
concentration at eroding sites in comparison to depositional sites. 

Soil erosion drastically influences not only lateral SOC 
distribution within a landscape but also vertical CO2 fluxes into 
the atmosphere [7,10]. Van Oost et al. [5] summarized at least 
three key mechanisms controlling the net flux of C between the 
soil and atmosphere: 1) dynamic replacement of SOC at the 
eroding sites [6]; 2) deep burial of SOC rich topsoils at 
depositional sites [4,13]; 3) enhanced decomposition of SOC 
because of the chemical or physical breakdown of soil during 
detachment and transport [7]. Particularly, the second and the 
third mechanisms should be susceptible to changes in the 
precipitation regime. 

A key uncertainty of erosion-induced C loss is C mineralization 
resulting from the breakdown of soil aggregates as a direct 
response to extreme precipitation [7,14,15]. During a given 
erosion event, rainfall leads to breakdown of aggregates and 
releases the encapsulated C due to flow shear and raindrop impact 
[14]. Some studies suggest that aggregates breakdown by raindrop 
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impact and wetting is mainly caused by initial fast slaking [16] or 
welding [17]. However, the extent of additional CO2 fluxes from 
breakdown of aggregates due to erosion is still largely unknown. 
Franzluebbers [18] estimated a 10-60% increase in CO2 evolution 
from various soils after breakdown of aggregates during 0-3 days. 
Polyakov and Lai [14] suggested that mainly the breakup of initial 
soil aggregates by erosive forces is responsible for increased CO2 
emission. However, conducting a set of rainfall simulation 
experiments, Bremenfeld et al. [19] recently suggested that interill 
erosion and associated soil aggregates breakdown have no 
prominent effect on soil respiration in situ. Therefore, effects of 
erosion-induced breakdown of aggregates on COg evolution need 
to be further assessed. 

Estimates of soil and SOC redistribution and associated GO2 
emissions show a large spatial and temporal variability. As field 
SOC and CO2 fluxes of soils under erosion strongly depend on 
temporsd variability of environmental conditions (e.g. location, soil 
management, initial soil moisture, and rainfall event characteris- 
tics) rainfall simulations under controlled laboratory conditions 
may help to shed light on C flux processes. Several rainfall 
simulation experiments have attempted to investigate soil erosion 
and associated SOC dynamics [20-24]. Jacinthe et al. [24] 
determined mineralization of SOC in runoff under no-tiU, chisel 
till and moldboard plow conditions with rainfall simulation 
approach. Van Hemelryck et al. [23] experimentally simulated 
three t)'pii:al agTi[:ulture erosion events to quantify CO2 emission. 
So far, however, there is no direct process assessment on 
combining effects of erosion, transport and subsequent deposition 
on C redistribution including vertical CO2 fluxes. Changes in 
SOC pools indicative for important mechanisms of SOC 
redistribution and differing in their stability against microbial 
decay are not well known. 

To get a better process understanding of soil erosion, transport 
and deposition on the redistribution and mineralization of SOC, 
the main objective of the present study was to determine SOC 
redistribution and a complete C budget of a loess soil affected by 
water erosion using a pseudo-rephcated rainfall simulation 
experiment under standardized conditions. The following pro- 
cesses were studied and considered in our C budget: 

(i) We determined SOC mineralization by measuring CO2 
emissions at different slope positions. 

(ii) We analyzed soil and C redistribution along the slope 
including potential export into aquatic ecosystems. We 
measured C enrichment in the redistributed sediment. In 
order to test the hypothesis that POC is preferentially 
eroded and exported into aquatic ecosystem we fractionated 
SOC by density into particulate organic C (free POC, C not 
bound to minerals) and mineral associated organic C 
(MOC). 

(iii) Finally, we analyzed concentrations of DOC in soil solutions 
at different positions of the slope and in runoff and 
determined above and belowground lateral DOC fluxes. 

Materials and Methods 

Ethics Statement 

The experimental station 'Proefboerderij Wijnandsrade' (The 
Netherlands) permitted access to their land and allowed for taking 
soil sample material from their cereal fields for the research carried 
out. 



Site Description and Sampling 

The loess soil was collected from an agricultural field with 
winter wheat in South Limburg (50°53'58. 42"N, 5°53'16. 23"E), 
The Netherlands in May 20 1 1 . South Limburg is part of the 
European loess belt and has a temperate maritime climate. This 
region has a mean annual precipitation of 825-850 mm [25] and 
a mean annual temperature of 10.2°C. The sampled soil has a sUty 
loam texture, and is classified as a Haplic Luvisol [26]. In the 
present study, the top 10 cm of the Ap horizon was collected and 
sieved over an 8 mm mesh to homogenize the soils and to keep 
aggregates intact as much as possible. Agricultural management at 
the sampled site is characterized by a potato-winter wheat-beet- 
winter wheat rotation. Soils are plowed 30 cm by a cultivator in 
spring and conventional tillage was applied in winter (including 
30 cm plowing). The basic physical and chemical properties of the 
used soil are shown in Table 1. 

Soil Analysis 

Field bulk density was estimated from undisturbed 100 cm ^ 
cores that were oven-dried at 105°C for 24 hours [27]. Grain size 
distribution of soils was obtained using a particle size analyser 
(Micromeritics, SediGraph 5100, Norcross, USA). Soil pH (1:2.5 
in H2O) was measured with a multi-parameter analyser (CON- 
SORT C832, Abcoude, The Netherlands). Soil water content was 
continually determined by a multi-channel Metallic TDR cable 
tester system [28]. Carbon and nitrogen (N) contents in bulk soils, 
sediments and density fractions were determined using a C and N 
analyser (Elementar VarioEL, Hanau, Germany). 

Experimental Design 

The erosion experiment was carried out using a 
1.25 mx3.75 m experimental stainless steel flume (Figure 1). 
The upper 1.75 m had a slope of 15° (upslope position) and the 
lower 2 m had a slope of 2° (downslope position). To assess the 
effects of erosion, transport and sufjsequent deposition on 
redistribution of soils and C along the erosion slope, the 
experimental flume was divided into three zones according to 
the positions of the slope and observed results of sediments 
redistribution (Figure 1): 1) the eroding zone, at the upper half of 
the upslope position; 2) the transport zone, at the lower half of the 
upslope position and the upper half of the downslope position; 3) 
the depositional zone, at the lower half of the downslope position 
of the flume. We used a static definition of the different zones as 
dynamic measurement locations would have disrupted the soil 
surface. We recognize that these zones can change during the 
event and between events and that during events in every zone 
also local deposition and re-entrainment will occur. 

The entire flume was subdivided into three parallel replicates of 
40 ±2 cm wide. The soil was laid on top of a 2 cm thick layer of 
inert quartz sand to allow water to drain away. On top there was a 
20 cm layer of soil on the upper (erosion) section where soil \vas 
supposed to erode and a thinner (5 cm) soil layer on the lower 
deposition section to ciUow for material deposition. On the 
transport section there was a gradual transition from 20 to 5 cm 
soil layer. While placing tfie air-dried soil it was compacted for 
every 2 cm, using a hammer and wooden piece of board 
(30x30 cm) to distribute the applied force. The compaction was 
such that it approached bulk density under field conditions (1.28 g 
cm ^). In addition, there were three controls. Three control 
buckets (diameter 34 cm) were filled with a 20 cm loess soil layer 
on top of a 2 cm quartz sand layer, similar to the main flume. 
These controf buckets were also compacted to the same bulk 
density. The buckets were placed next to the flume so that they 
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received the same rainfall as well, but no lateral displacement of 
soil material took place. 

The soU layer was pre-wetted to an initial standard moisture 
contents (Table 2) in 10-15 min to initiate runoff generation prior 
to commencing the real rainfall experiment. Four 18-minutes 
rainfall events were carried out at a monthly time interval. 
Measurements were carried out every 2 minutes during rainfall 
simulation. Rainfall was simulated with two nozzles (Lechler 460 
788) applying at 1600 hPa demineralized water using an average 
rainfall intensity of 41. 8± 1.9 mm h~ . A rainfall event with this 
intensity and duration of 18 minutes has a return period of about 2 
years [29]. Mean drop size of the applied rainfall was 2.0 mm 
(D5o = 2.0). With an average falling height of 1.8 m, the kinetic 
energy applied on the soil surface was 12. 5J m~^ mm~'. 
Demineralized water was used instead of tap water to prevent 
flocculation problems with dispersible soil material [30,31]. As the 
total load of ions in rainwater is very low (the annual average 
electrical conductivity EC25 was below 20 \lS cm ' at the official 
Dutch sampling site Beek [32], about 10 km from the soil 
sampling site) the physico-chemical impact of demineralized water 
on soil particles is considered to be the same as for rain water. The 
temperature was kept as constant as possible (18. 1 ±0.9°C). 

Sampling during Erosion Experiments 

Sediment traps were installed in the middle of the eroding, 
transport and depositional zones respectively with entrance of the 
traps at the upslope side and at the same level as the soil surface to 
capture mobilized sediment in overland flow (Figure 1). The 
sediment traps were modified 12 ml Polypropylene screw cap 
tubes (Greiner Bio-One GmbH, Frickenhausen, Germany). The 
traps had a small diameter to minimize disturbance to overland 
flow and resulting erosion patterns. An opening in the side was 
made to collect mobilized sediments. The sediment traps were 
sampled every two minutes and the collected materials were 
transferred to containers, oven-dried at 35°C, weighed and later 
analysed for C and N contents. 

Runoff and sediments were collected from weirs at the end of 
the flumes at 2-min intervals once continuous runoff had 
developed. Total runoff was collected using a polystyrene gutter 
that was installed at the lowest part of the experimental flume. The 
contents of the flume were then pumped into V-notched botties to 
measure flow rates using a simple siphon pump made of Tygon R- 
3603 tubes (Saint-Gobin, Courbevoie, France). The lower end was 
constrained to 4 mm diameter to provide a constant flow velocity, 
without risking clogging by larger soil particles and keeping effects 
on the aggregation of the sediments limited. The V-notched 
bottles overflowed into sampling boxes which were replaced every 
two minutes or when the sampling box was fuU. 

At the lowest end of the flume, three holes per replicate flume 
were present at the level of the sand drainage layer to collect 
through flow. Through flow was defined as the lateral under- 
ground flow in contrast to the overland flow. The holes were 
covered from the inside by a 63-|J,m stainless steel mesh allowing 
water to pass through, but to prevent clogging up. On the outside 
of the walls attached tubes drained into botties, similar to the 
runoff setup. 

Sampling after Erosion Experiments 

Density Fractionation of Bulk Soils and Exported 
Sediments. After four rainfall events the 0—2 mm topsoUs at 
the eroding, transport and depositional zones and sediments 
exported during the first and fourth events were fractionated into 
three fractions by a sodium polytungstate (NaPT) solution with a 
density of 1.6 g cm the free light fraction (fLF) which consisted 
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Figure 1. Photographs of the experimental setup and sampling locations along the experimental flume. It included the eroding, 
transport and depositional zones of the flume. A shows the lateral view; B shows the vertical view. 
doi:1 0.1 371 /journal.pone.0096299.g001 



of large, undecomposed or partly decomposed root and plant 
fragments, the light fraction occluded in aggregates (oLF) and the 
heavy fraction (HF), which was associated with minerals [33,34]. 
Soil organic C in fLF, oLF and HF are defined as fPOC, oPOC 
and MOC, respectively. Particulate organic C (POC) is the C not 
bound to soil minerals including both fPOC and oPOC. The 
oPOC represents C sequestered in aggregates. Methods and 
procedures were followed as described in Cerli et al. [34]. AU 
fractions were freeze-dried, homogenized and later analysed for C 
and N contents. Density fractionation was done in triplicate. 

Dissolved Organic Carbon (DOC) 

To investigate dynamics of dissolved organic C in different soil 
depths and positions as affected by soil redistribution, soil moisture 
samplers (MACRO RHIZON 19.21.35, 9 cm porous, 4.5 mm 
OD, 0.2 |J.m, Wageningen, The Netherlands) were inserted in the 
eroding, transport and depositional zones of the flume. Each 
sampler was connected to a syringe (50 mL) to collect the soil 
solution. At the eroding and transport zones of the flume, soil 
solutions were collected at 4 cm and 9 cm depths. At the 
depositional zones soil solutions were sampled at 4 cm only 
because of the thinner soil layer on the lower deposition section. 
Soil solutions were sampled twice per week during the first week 
immediately after one rainfall event because of higher soil water 
moisture. As the soil dried, soil solutions were collected once per 
week. Concentrations of DOC were determined by a TOC 
analyser (TOC-V CPH, Shimadzu, Kyoto, Japan). 

Soil CO2 Efflux Measurements 

Soil respiration was measured using a Portable Gas Exchange 
and Fluorescence System (LT 6400XT; LICOR Biosciences, 
Lincoln, NE USA). In order to enhance the comparability of data, 
most CO2 efflux measurements were conducted in the afternoon 
between 17:00 and 19:00 at local time in PVC collars (10.2 cm in 
diameter and 7 cm in height). Soil CO2 efflux was determined 
before and after each rainfall simulation event. As the 7 cm high 
collars, necessary for the CO2 efflux measurements, would strongly 
affect the overland flow and erosion patterns during the rainfall 
event, the 7 cm high collars were replaced by smaller collars (same 
diameter but 1.5 cm tail). These were inserted at exactly the same 
place, to temporary fiU the imprint of the high coUar in the soil 
surface. The top of the coUar was placed exactly equal to the soil 
surface, to minimize the disturbance of the sampling location by 
the CO2 measurements but still enabling to measure the CO2 
efflux exactly at the same position later on. Overland flow was 
possible and erosion, transport and deposition processes at the 
surface of the area used for measuring CO2 were hardly affected 



by this strategy. Two to three measurements per site (i.e. per coUar) 
were carried out each time. The number of replicated measure- 
ments per collar depended on the variation after the first two 
analysis with an additional measurement if the relative deviation of 
the second one was larger than 10%. Additionally, pre- 
experiments were carried out using the same loess soil to test 
impacts of soil depth on soil CO2 efflux. In these experiments the 
CO2 efflux was measured in columns with increasing soil thickness 
under constant soil moisture and temperature conditions. Results 
showed that soil depth did not have significant effect on soil 
respiration per soil weight up to a depth of 30 cm (data not 
published). Based on these results, aU data measured in different 
experimental zones, and control soils, having different soil depths, 
were corrected to 20 cm soil layers in order to directly compare 
effects of erosion, transport and deposition on CO2 effluxes. 

Erosion-induced Carbon Budget 

Fluxes of SOC and DOC were calculated by multiplying 
concentrations of SOC and DOC with the volume of the overland 
flow. Other parameters were calculated as follows: 



Carbon enrichment ratio (CER)- 



C, 



C control soils 

Total C losses = Lateral C exported 
+ Vertical CO2 emission 

Lateral C exp orted = SOC exported in overland flow 
+ DOC exported in overland flow 
+ DOC exported in through flow 

Net additional CO2 emission 

= CO2 emission from soil in flume 

— CO2 emission from soil in control treatment 



Based on the 4-month data we calculated annual C fluxes by 
linear extrapolation making comparisons with the literature easier. 
However, the shortcomings of such budgets based on short-term 
laboratory experiments only are obvious. 

The definition of C source and sink areas for calculating the C 
budget was based on two experimental observations. After the 
fourth rainfall event, soil layers with relocated materials were 
clearly visible in the flume, particularly in the downslope part of 
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the depositional zone (Figure IB). In addition, we found that SOC 
was significantly depleted in the transport zone comparing with 
controls soils (cf. section results). Based on these two observations, 
the eroding and transport zones were defined as the C source area 
and the depositional zone and the runofi' leaving the flume 
(exported into aquatic system) were defined as the C sink area. We 
calculated an erosion-induced SOC budget for the four rainfall 
events over the entire period using a mass balance approach (i.e. 
source = sink area). Changes in C distribution between the 
density fractions were appropriately considered by using the data 
of the original soil for the source area. This approach enabled us to 
include any changes in C redistribution between density fractions 
induced by erosion. 

Statistical Analyses 

Differences in C enrichment ratios, amounts of sediment 
exported and DOC concentrations in overland flow were tested 
with one-way ANOVA and the Post-hoc Duncan test to 
diSerentiate between individual differences. The difference of 
CO2 effluxes measured in the 4-month period at eroding, 
transport and depositional zones of the gutter was tested by 
repeated measurement ANOVA. Averaged CO2 efflux in different 
experimental zones was compared using a one-way ANOVA. For 
all tests, a significance level of -P= 0.05 was set using the Post-hoc 
Duncan test, unless otherwise indicated. The relationship between 
cumulative CO2 emission and DOC concentration was tested by 
two-tailed Pearson test. All statistical tests were performed using 
SAS software (Version 8.1) and SPSS (IBM Statistics 20). 

Results 

Loss of Sediment and Carbon Enrichment Ratios in 
Overland Flow 

Total sediment losses in the overland flow increased during the 
course of the experiment from 9.5 g m ^ in the first event to 
3 1 .0 g m ^ in the fourth event (Figure 2). During the first rainfall 
event the average sediment concentration was 1.1±0.2 g L~' and 
doubled to 2.3±0.8 g L 'in the fourth event. 

Carbon enrichment ratios (CER) of sediment loads of overland 
flow trapped at the eroding, transport and depositional zones of 
the flume ranged from 0.8 to 2.9. The CER was significandy 
higher at the depositional zone compared to those of the eroding 
and transport zones (Table 3). Carbon enrichment was even 
stronger in the sediments of the runoff with CER between 1.3 and 
4.0. Carbon enrichment ratios decreased with increasing concen- 
trations of suspended sohds in the overland flow (Figure 3). 
Concentrations of suspended solids were smaller at the beginning 
of each rainfall event, resulting in larger C enrichment but also in 
larger variation of the data. 

Preferential Erosion and Deposition of Organic Carbon at 
the Soil Surface 

After four rainfall events, a thin sedimentation layer was present 
in the depositional zone (approximately 2 mm thick) without any 
layering. However the depositional zone clearly showed patterns of 
deposition of finer grained materials along the flow lines of 
overland flow and the whole lower part of the gutter. Soil organic 
C concentration (mg ' g soU) of the surface soil (2 mm) decreased 
by 6.0% at the eroding zone and increased by 3.9% at the 
depositional zone if compared to control soils (Table 3). Never- 
theless, soil organic C concentration did not differ significandy 
between the control, the eroding, transport and depositional zones 
of the gutter. Also the relative distribution of C in density fractions 
of the soil was not affected by soil erosion. Most of the C (86% to 
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Figure 2. Average total eroded sediment per rainfall event exported by overland flow and carbon enrichment ratios (CER) during 
four rainfall events. 
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91%) was found in the heavy fraction, i.e. mineral associated 
organic C (MOC; Table 3). The rest was almost equally 
distributed between the free light fraction (particulate organic C 
in free light fraction = fPOG) and the fraction occluded within 
aggregates (oPOC). The free light fraction was significantiy 
enriched in C at the transport and the depositional zone whereas 
the occluded light fraction (oLF) was depleted in C at the eroding 
zone (Table 3). The heavy fractions of surface soils in the flume did 
not significantly change in C contents. 



In the sediments, all fractions were strongly enriched in C with 
the largest enrichment in the free light fraction. This C enrichment 
was smaller in the occluded and smallest in the heavy fraction and 
also decreased from the first to the last event (Table 3). However, 
the C content of the heavy fraction of the sediments (first event) 
was more than double the C content of the heavy fraction of tlie 
control soil (Table 3). 



4.5 




1.0 -I 1 1 1 1 1 

0 1 2 3 4 5 6 

Suspended solid concentration (g kg"^) 

Figure 3. Relationship between carbon enrichment ratio (CER) and suspended solid concentration (SSC) in the overland flow. 
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Soil CO2 Efflux 

All measured CO2 efflux rates for the whole experiment ranged 
from 0.12 to 4.34 g C m~^ day"' (Figure 4). During the entire 
experimental period, rates of CO2 emissions exhibited a similar 
behaviour in the eroding, transport and depositional zones and the 
non-eroded control with a sharp initial increase immediately after 
each rainfall event, followed by continuously decreasing rates 
thereafter. Rates of CO2 efflux significantly decreased with time 
during the four events (P= 0.001). The spatial and temporal 
variability of CO2 efflux rates was larger in the first rainfall event 
than during the other events. 

The largest mean CO2 efflux was observed in the transport zone 
during the first three rainfall events (Figure 5). In the fourth event, 
■S ^ however, the depositional zone had the largest mean CO2 efflux. 

S g The relative differences of the mean CO2 efflux between the 

0 ^ depositional and the eroding zones increased during the course of 
T c the whole experiment and became significant in the fourth event, 
u 5 Cumulative CO2 fluxes in the eroding, transport and deposi- 
l" a tional zones ranged from 80 to 180, 1 16 to 317, and 146 to 204 g 

£ C m~^ yr~'; respectively. The largest mean cumulative CO2 

g_ ^ fluxes (221 g C m ^ yr ') were observed in the transport zone. 

^ Mean cumulative CO2 fluxes in the depositional zone (181 g G 

75 2 m~^ yr~') were significandy larger than those in the control soils 

c t (P=0.02) while CO2 fluxes in the eroding zone were similar in 

— "2 comparison to the control. The total losses of C as CO2 emission 

.9 g during the entire experiment accumulated to 1.8 to 2.9% of total 

J '-^ soil organic C stocks. 

1 I 

^ Y Dissolved Organic Carbon (DOC) 

^ 5 Concentrations of DOC in soil solutions at eroding, transport 

^ S and depositional zones ranged from 7.1 to 25.9 mg L~ during 

J s four rainfall events (Figure 6). In the shallow soil (4 cm depth), 

° _g mean concentration of DOC decreased in the following order: 

.S" g transport zone (15.1 mg L ') > control soils (14.3 mg L ') > 

_^ 'g depositional zone (12.3 mg L~') > eroding zone (11.8 mg L~'). 

§ S However, only DOC concentrations in the depositional and 

_g g eroding zones were significandy lower than those in the transport 

o y I zone and the control. Mean concentrations of DOC in the deeper 

g § >g soil ( 1 0 cm) were almost equal as in the shallow soil and decreased 

S o ^ in the following order (not statistically significant): control soils 

^ S I (16.8 mg L ') > transport zone (15.2 mg L ') > eroding zones 

uT^^ (12.3 mg L '). Concentrations of DOC in soil solutions of both 

S c -o depths showed distinct temporal patterns in all zones of the gutter, 

c tj ™ They increased at the beginning of each rainfall event, then 
decreased and increased again with time. This trend was less 
obvious during the first rainfall, 

y CTi S- Concentration of DOC in overland flow remained constant 

Q. E 



o = -2 



w- 9 .! 



o ^ 'ft S during each single event, ranging from 0.3 to 8.3 mg L ' and 
significandy decreased from the first to the third rainfall event 

f ^ 5 o (means of the four rainfall events: 7.2±0.4, 2.6±0.4, 0.9±0.7, 

g, j^ >-:5 3 0. 7 ±0. 4 mg L~', no further data shown). Cumulated DOC fluxes 

° oi c o c transported by overland flow were on average 0.23 g C m ^ yr ' 

3 "i^ ^ ? u o (Figure 7). The amount of C exported as DOC by overland flow 

S .i "i ? 5 was small, accounting for 0.014% of the total SOC stocks in the 

^ liSuS flume. Fluxes of DOC in through flow (i.e. 0.002% of total SOC 

^ 'S 2' S § stocks) were significantly smaller than overland flow. 

N c c 2 'ft ^ 

.2 3 2 s £ 2 .2^ Discussion 

£ Si § S I §.=5- Preferential Transport and Deposition of Organic Carbon 

"^■^cccj^ As expected from the literature [8,12], the soil of the eroding 

^ ° -9 -2 -S ™ 5 zone was depleted in C whereas the soil of the depositional zone 

— ^ ^ ^ '-J '^ ■q and the sediments of the overland flow were enriched in C after 



the four rainfall events (Table 3). The results of the density 
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Figure 4. CO2 efflux at different zones of the gutter and control soil during four rainfall events. Solid line + circle represents the eroding 
zone; dotted line + circle represents the transport zone; solid line + triangle represents the depositional zones; and dotted line + triangle represents 
control soils. Values are mean± standard error of three replicates. 
doi:1 0.1 371 /journal.pone.0096299.g004 

fractionation clearly showed a large loss of C occluded in 
aggregates in the eroding zone, which was accompanied by an 
enrichment of C in the fPOC fraction in the other zones of the 
flume and in overland flow (Table 3). We assume that the 
disruption of macro-aggregates by raindrop peeling [35] and 
aggregate welding and development of a structural crust [17] 
resulted in the liberation of fPOC, which was preferentially 



transported [8] . The disruption of macro-aggregates will result in 
the release of micro-aggregates (smaller than 250 |J.m) from the 
macro-aggregates too. The C content of micro-aggregates within 
macro-aggregates is usually larger than that of macro-aggregates 
[36-38]. The release of such small aggregates and selective 
transport of small aggregates with low density [39] could be the 
reasons for the observed significant C enrichment of oPOC in 
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Figure 5. Mean cumulative CO2 emission at the eroding, transport and depositional zones and control soil. Different capital letters 
mean significant difference at a single rainfall event between the different zones. Values are mean± standard error of three replicates. 
doi:1 0.1 371/journal.pone.0096299.g005 
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Figure 6. Dissolved organic carbon (DOC) concentrations. DOC solutions were collected at 0-4 cm and 0-9 cm depths of the eroding, 
transport and depositional zones of the flume during four rainfall events. Solid line + circle represents the eroding zone; dotted line + circle 
represents the transport zone; solid line + triangle represents the depositional zones; and dotted line + triangle represents control soils. Values are 
mean± standard error. 
doi:1 0.1 371 /journal.pone.0096299.g006 



sediments ranging from 2.3 to 3.2 (Table 3). However, we did not 
study aggregate stability and the detailed processes resulting in 
breakdown of the aggregates neither the related preferential 
erosion, transport and deposition of different sizes of aggregates 
and particles. That should be done in follow-up experiments. 

The calculated mass balance of the experiment illustrates the 
disruption of aggregates in the eroding zone and the redistribution 
of C from aggregates to fPOC with an erosion-induced 
accumulation of fPOC in the sink area of 0.24 g C (Table 4). 
This accumulation is equal to an increase in fPOC by 48% 
comparing the source and the sink area. One logical source of this 
additional POC would be C occluded within aggregates in the 
eroding zone at the beginning of the experiment. This large 



accumulation of fPOC in the sink area contributed to the observed 
relative increase by 6% (Table 3) in the SOC content of the first 
two mm layer of the depositional zone. 

Particulate organic C already present in the soils and formed by 
disruption of aggregates (cf above) was preferentially eroded and 
transported by overland flow as indicated by the largest CER ratio 
of the density fractions in any of the sampled soils and sediments. 
Per definition, fPOC is the lightest fraction, not associated with 
minerals and therefore easier to be translocated by water than soil 
particles with a higher density [23,40,41]. The high C enrichment 
of mineral-associated organic C (MOC) in the sediments of the 
overland flow suggested that water erosion separated the whole 
soil particles according to their density (Table 3). This fraction- 
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Figure 7. Conceptual diagram illustrating the total carbon budget as affected by soil erosion, transport and deposition in the four 
months rainfall simulation experiment. Fluxes were calculated on an annual base (interpolated from the 4-months experiments). The values 
were expressed as mean values and standard error of three replicates. 
doi:1 0.1 371 /journal.pone.0096299.g007 



ation occurred between the different density fractions. Increasing 
C concentrations (MOC<oPOC<fPOC; Table 3) resulted in 
increasing CER of the sediments in the same order. 

A significant portion of the eroded and transported C enriched 
sediment was not retained in the downslope parts of the 
depositional zone and was exported by overland flow and left 
the flume (Table 4). Particularly the weakly decomposed C of the 
fPOC should be a readily available C and nutrient source for 
aquatic organisms [16,42] contributing to CO2 emission from 
aquatic ecosystems. This process linking terrestrial and aquatic 
systems cannot be neglected for modeling the C cycle and has to 
be studied in more detail in future. 

Relationship between Erosion Rate and Carbon 
Enrichment 

The inverse, non-linear relationship between the erosion rate 
and C enrichment of the sediments we found (Figure 3) is in 
agreement with previous studies [21,40,41,43]. This inverse 
relationship is the result of increasing sediment concentration in 
the overland flow during each single event and from the first to the 
fourth rainfall event. One of the most important reasons for this 
relationship should be the breakdown of macro-aggregates by the 
raindrops as already discussed to be the main reason for the 
preferential erosion of fPOC [8,10,35]. This process should be 
particularly important at the beginning of each rainfall event 
because rewetting of dry soils results in the disruption of aggregates 
and the release of organic matter [44]. It is also reasonable to 
assume that the importance of this process wiU decrease with 
increasing number of rainfall events. Heavy rainfall causes 
compaction, welding and crust formation resulting in reduced 
infiltration and increased erosion and suspended solid concentra- 
tion with time [15]. The preferential removal of C enriched soil 
wiU result in C enriched sediments particularly at the beginning of 
the experiment where the erosion rate was stiU small. After 
removal of this soil enriched in C, the erosion rate increases 
because of decreasing infiltration and generation of more overland 
flow. That will result in even increasing erosion rates because soils 
are less protected by organic matter and aggregation. In field 
situations, the relationship between erosion rate and C enrichment 
might be weaker because of continuous above and belowground C 
input and its positive effect on aggregation [45]. 



Decreasing C enrichment with large erosion rates, i.e. 
increasing sediment concentration, indicated that an increasing 
erosion rate does not result in proportionally increasing C losses. 
However, this does not mean that more severe erosion events lead 
to less impact on soil C. Very strong erosion events will translocate 
large amounts of C. However, this C might be better protected 
against further mineralization after deposition because G is mostiy 
deposited as mineral associated C (Table 4). The C loading of 
mineral surfaces should be low as well, resulting in a more efficient 
stabihzation against microbial decay [46,47]. In addition, long- 
term erosion-induced C sequestration or depletion might depend 
on the precipitation frequency and intensity. 

Soil CO2 Effluxes 

This study provides new data and insight on C decomposition 
under controlled conditions in an artificial landscape setting at 
eroding, transport and depositional positions allowing for a better 
process understanding. Although we did not scale up our results to 
the landscape level, it is important to know whether the fluxes 
measured do compare with observed field measurements and 
make any sense, also in comparison with previous indirect 
measurements on eroded sediments and soU profile investigations 
[5,12]. In the present study, CO2 efflux rates measured (0.12 to 
4.34 g C m ^ day ') were in the range of soil respiration rates 
from agricultural loess soils [5,19]. Initial increases of CO2 
emissions immediately after each of our rainfall events might be 
explained by the increase in microbiological activity after re- 
wetting the dry soil and/or increased bioavailability due to 
aggregate breakdown [44] . Aggregate breakdown and subsequent 
exposure of previously encapsulated SOC provide substrates for 
microbial decomposition [8]. The re-wetting effect was particu- 
larly important after the first event and decreased during the 
course of the experiment. This is in line with a decreasing capacity 
of soil to release C from aggregates over time [48,49]. 

Transport of topsoU and associated C influenced SOC 
decomposition rates at the different positions of the artificial 
slope. The small cumulative CO2 emission from the eroding zone 
should be the result of the observed preferential removal of C 
enriched materials (i.e. higher CER of the sediments at the 
depositional zone in comparison to the eroding zone, Table 3), 
which was either preferentially deposited or left the gutter. This 
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preferential removal of more labile C (POC: fPOC and oPOC) left 
behind less C, which was relatively more stable (Table 3). In turn, 
accumulation of labile C fractions in the depositional zone 
contributed to an increasing difference in cumulative GO2 
emission between the eroding and depositional zone over time 
(Figure 5). 

The cumulative CO2 emission was significantly and positively 
related to fPOC (R^ = 0.94; P= 0.03) illustrating the more labile 
character of this SOC fraction. The large accumulation of fPOC 
probably explained that the transport zone had the largest 
cumulative CO2 emissions (Table 3 and Figure 7). Larger CO2 
emissions from the depositional zone were anticipated because the 
deposited labile C (fPOC, oPOC) could be used as substrate and a 
source of energy for microbial respiration [50,51]. Although CO2 
emissions were large at the depositional zone, the SOC content 
increased by 6% in comparison to the control soils after four 
erosion events. Obviously, parts of the eroded and deposited SOC 
were preserved (Figure 5). 

Considering all positions of the slope, mean DOC concentration 
in near-surface layers was positively correlated to median soil CO2 
efflux rate {P= 0.02). The largest CO2 efflux was accompanied by 
largest DOC concentration in the transport zone - a second 
parameter (first fPOC) explaining the large CO2 eflQux in this 
zone. Creed et al. [52] found that substrates (i.e. DOC) in the 
near-surface soil were strongly related to median soil CO2 efflux. 
Considering each position of the slope separately, median soil CO2 
efflux rates were not significantly related to mean DOC 
concentrations at the eroding (P=0.18), transport (-P=0.49) and 
depositional zones {P— 0.22). However, DOC was significantly 
correlated to the median soil CO2 efflux rate in the control soil 
{P= 0.05), which indicated DOC could be mineralized during the 
experimental period. Thus, DOC dynamics could not explain the 
observed additional C decomposition at the depositional zone. 
This might suggest a fast turnover of DOC or/ and a direct use of 
POC by the microbial community. 

Total Carbon Budget 

We estimated an erosion-induced C loss of 53 g C m~^ yr~' 
calculated as the sum of net erosion-induced CO2 emission, C 
losses by overland flow (C in sediments, DOC) and through flow 
(cf materials and methods, figure 7). During the entire experi- 
mental period, the averaged SOC fluxes leaving the flume with 
overland flow were 18 times larger than DOC fluxes including 
lateral fluxes by the through flow (Figure 7). Fluxes of DOC 
(0.26 g C m ^ yr ') were rather low particularly due to decreasing 
DOC concentration during the experiment, i.e. with increasing 
number of events. Fluxes of sediment associated C were equivalent 
to 8.9% of the erosion-induced C loss while DOC fluxes were 
equivalent to 0.5% of those C losses. Therefore, sediment 
associated C played a much larger role than DOC in the 
erosion-induced linking of terrestrial and aquatic ecosystems. 
However, erosion-induced DOC flux should not be neglected 
because DOC might be particularly important for aquatic food 
webs [42,53]. 

Erosion-induced CO2 emission was the dominant form of C 
loss, representing 90.5% of erosion-induced C loss. Based on the 
assumption made (cf material and methods), 1.5% of total C 
redistributed (deposited C at the depositional zone plus C exported 
to aquatic ecosystems) was mineralized to CO2 (Table 4 and 
Figure 7). Previous estimates of decomposition of eroded SOC 
showed large variations, ranging from 0 to 100% (e.g. 
[7,14,15,23]). In modelling studies the assumption is often used 
that at least 20% of the eroded SOC is decomposed as a 
consequence of soil erosion [7,54]. Our measured values were 
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much smaller than this conventional view of erosion effects on the 
C cycle. Palyakov and Lai [14] estimated 8% of SOC displaced by 
erosion had a potential to be mineralized and van Hemelryck et al. 
[23] estimated mineralization of 2% to 12% of the eroded SOC in 
a loess soils using laboratory rainfall simulation experiments. We 
propose three main reasons for the large difference. Firstly, the 
effects of the disruption of aggregates on extra CO2 efflux were 
relatively short-lived [48] . Secondly, C stabilization as affected by 
soil erosion and deposition might be underestimated in the 
previous assumption [5,8]. Thirdly, the artificial slope was 
relatively short in our experimental setting in comparison to the 
field, which may result in an underestimation of transport effects 
on C mineralization. 

Conclusions 

Rainfall simulation experiments are a useful approach to 
determine the role of soil erosion for the C cycle. The data of 
our 4-months experiment were comparable to field situations, 
despite of well-known shortcomings of laboratory approaches. The 
erosion rate was estimated to be 2.1 mm yr (26 t ha yr ), 
which was comparable with estimations in this region based on 
field data [55] . Also C enrichment of exported sediments and soil 
CO2 efflux were in the range of field measurements. 

Erosion-induced CO2 emission was the dominant form of C 
loss, representing about 90.5% of the erosion-induced C loss. In 
addition, a considerable amount of C rich sediments (265 g m ^ 
yr was laterally exported by overland flow. Carbon associated 
with sediments was the main form of erosion-induced lateral C loss 
and not DOC. This exported C plays an important role in the 
connection of terrestrial and aquatic: ecosystems. 

In our experiment, this redistribution of C rich materials 
resulted in a net additional CO2 emission during transport and 
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